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In hopping magnetoresistance of doped insulators, an applied magnetic field shrinks the electron 
(hole) s-wave function of a donor or an acceptor and this reduces the overlap between hopping 
sites resulting in the positive magnetoresistance quadratic in a weak magnetic field, B. We extend 
the theory of hopping magnetoresistance to states with nonzero orbital momenta. Different from 
s-states, a weak magnetic field expands the electron (hole) wave functions with positive magnetic 
quantum numbers, m > 0, and shrinks the states with negative m in a wide region outside the 
point defect. This together with a magnetic-field dependence of injection/ionization rates results in 
a negative weak-field magnetoresistance, which is linear in B when the orbital degeneracy is lifted. 
The theory provides a possible explanation of a large low-field magnetoresistance in disordered 
TT-conjugated organic materials (OMAR). 
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As is well known the magnetoresistance (MR) in the 
hopping regime is caused by a strong magnetic field de- 
pendence of the exponential asymptotic of the bound 
state wave function at a remote distance from a donor (or 
an acceptor) [l|. In the case of the Coulomb potential, 
when the wave function is spherically symmetric in the 
absence of the magnetic field, it becomes cigar-shaped 
squeezed in the transverse direction to the field [l|, 
This leads to a significant decrease in the overlap of the 
wave-function tails of two neighboring donors, and hence 
to a significant increase of resistivity (positive MR). An 
exponential positive MR in sufficiently strong magnetic 
fields is used to be a hallmark of the hopping conduction. 
On the other hand, there is anomalous (negative) MR 
observed in some hopping systems, for instance in amor- 
phous germanium and silicon. It has been attributed 
to magnetic-field dependence of spin-flip transitions be- 
tween sites when some fraction of them has a frozen spin 
Q, and/or to an increase of the density of localised states 
due to the Zeeman energy shift, fisB Q. This negative 
MR is used to be small (much less than 1%) even in rel- 
atively high magnetic fields of about 1 Tesla. 

In inorganic and organic insulators lattice defects such 
as vacancies, intcrstitials, excess atoms or ions and other 
"impurities" often localise carriers with a finite momen- 
tum rather than in the zero-momentum s-states. Here 
we extend the conventional theory of magnetoresistance 
[l[ to hopping via non-zero momentum orbitals. Quite 
remarkably this renders a giant weak-field magnetoresis- 
tance, which is negative. Moreover, if the orbital degen- 
eracy is lifted due to a broken time-reversal symmetry [sj 
with or without net magnetization, the negative MR is 
linear in B. 

The Schrodinger equation for the impurity-localised 
carrier wave function V'(r) can be written in the inte- 
gral form using the Green function, G{r,r'; E), (GF) of 



the Bloch electron in a magnetic field, 

i;{r) = - / drG(r, r'; E)V„np{r'Mr') (1) 



where E is the energy and Vimp(r) is the impurity po- 
tential. We consider first a two-dimensional (2D) system, 
such as a thin film in the magnetic field, perpendicular 
to the surface of the film. Generally GF is expressed as 
a sum over wave-functions of the 2D Bloch electron in a 
rational or irrational magnetic field with the Hofstadter's 
butterfiy eigenvalues @. In a weak magnetic field with 
hujc = heB/mt much smaller than the bandwidth, the 
effective band-mass (nib) approximation is sufficient, so 
that one can use the 2D free-electron GF in the magnetic 
field 
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where (j) and 0' are azimuth angles of ~^ and ~f^' respec- 
tively, / = {h/eBY^'^ is the magnetic length, T{a) is the 
gamma-function, and U{a, b, z) is the Tricomi's confluent 
hypergcomctric function well-behaved at infinity, z — >■ oo, 
for negative E [8]. Here a — 1/2 — {E ^ p,BB)/hu}c, 
where =p corresponds to spin up/down, respectively. Us- 
ing Eqs.dH [2]) one finds the wave function, ^p{r) — 
Fm(p) exp(im0), at ro ^ p ^ l^/ro as 

F^{p, h) oc pl™le-(«'')'''/«r(a);7 [a, I, [npfh/A] , (3) 

where rp is the radius of the impurity potential, k — 
{2mbeoy^^ /h is the inverse localisation length of the 
zero-field state with the ionisation energy eo, b = B/Bq 
is the reduced magnetic field with Bq — hK^/2e, and 
p is the distance from the impurity. While Fm{p) cc 
p^"^^G{p,0; E) is strictly applied to any finite-range 
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Vimpip), it works for the infinite-range Coufomb poten- 
tial at large distances as well Neglecting a small 
diamagnetic correction (quadratic in 6 <C 1) yields E — 
—eQ + hiOcm/2±fiBB where m = 0, ±1, ±2, ... is the mag- 
netic quantum number of the localised state, so that 



1 
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To elucidate the magnetic field dependence of the 
bound state we expand the solution, Eq. (|3]) , in powers of 
b making use of the integral representation of U (a, 5, z) 
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Replacing t with x — t/a and (1 + ax)"- 
(ax)"exp[l/x- l/(2aa;2) + l/{?,a^x^) - ...] yields 



(5) 



with 



&^[i^2M-i^oM] + 0(62), (6) 

where Kn{x) is the modified Bessel function. Eq.® is 
applicable in weak magnetic fields, 6^1 far enough but 
not too far from the point defect [up <!/&). In a wide 
region 1 <^ Kp < 1/h one can use the asymptotic [8] of 
Kn{x) w (7r/2x)i/2exp(-x)[l + (4^2 - l)/8z] to get a 
leading magnetic-field correction to the wave function, 
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Remarkably, the correction is linear in B for any non-zero 
m and could be large, if the bound state is sufficiently 
shallow and/or the hopping distance is large enough. It 
is positive for positive m and negative for negative m. 
The unusual expansion of the wave function with the 
positive m originates in the linear lowering of the ionisa- 
tion energy due to the orbital magnetic moment in weak 
magnetic fields. On the contrary, the states with nega- 
tive m shrink because their ionisation energy increases 
with the magnetic filed. As shown below the state ex- 
pansion/shrinking in the region Acp < l/b causes a linear 
MR in the weak magnetic field, which is negative or pos- 
itive depending on the particular orbitals involved in the 
hopping. 

In the case of a finite-range impurity potentials, Vij, 
the standard expression for the hopping integral 
Uj =< i \ Vj \ j > - < i \ j >< i \ Vj \ i >, yields 
tij (X Fjn{p), where p is the distance between two hop- 
ping sites i and j, which is assumed to be much larger 
than the localisation length. The hopping conductance 
is proportional to the hopping integral squared. Hence, 
the magnetoresistance to the hopping transport via par- 
ticular m-orbitals is found as 
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FIG. 1: (Color online) Hopping magnetoresistance in m = 2 
channel as a function of the reduced magnetic field B/Bo for 
a few hopping distances. 



where p = dxx^''"^^KQ{xy and q — 

dxx'^^~^^Fm{x,by accounts for some weak field 
dependence of the norm of the wave function. The same 
expression is also applied when hops take place between 
deep donor (acceptor) levels and shallow levels with the 
emission or absorption of phonons, as in the resistor 
network model of Ref . [§] . 

Here and below we take p as an average distance be- 
tween defects in the transverse direction to the field, so 
that p oc A^"^/^ in 2D systems, where N is the density of 
defects [l^l- FigH] represents hopping MR in the d-wave 
orbital state with m — 2 as a function of the reduced 
magnetic field for a few distances between defects. As 
one can see from the figure the anomalous expansion of 
the bound state results in the negative MR in weak mag- 
netic fields, b < [np)^^ followed by the positive MR in 
stronger fields. The absolute value of the negative MR 
could be in excess of 50 %. The negative (for m > 0) and 
positive (for m < 0) MRs are linear in B at relatively 
low B. 

Let us now extend the theory to a 3D system. To 
incorporate the free motion in the z direction one 
can replace E with E — pl/2mb and perform the 
Fourier transformation in Eq.®, G'3D(r, r'; £■) = 

{l/2^h) JZ^ dp.e^P^(-'-^)G2D (7^, 7^;; E - pl/2mt) . 
Replacing t in the integral representation of the confluent 
hypergeometric function, Eq.®, with t = l/(exp(x) — 1) 
and integrating over one obtains 
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Expanding the exponent in the square brackets in Eq. ([9]) 
up to the third power in x and performing the integration 
by the saddle-point method we finally obtain the asymp- 
totic of the 3D wave function, ^(r, &) cx pl'"lG3£)(r, 0; E) 
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FIG. 2: (Color online) Negative OMAR (%) at room temper- 
ature in ITO/PEDOT/AlqZ/Ca device at bias voltage 14 
V (symbols, Ref.Qll) described with Eq.((Tl| with r = 1.8, 
Bs = 64 mT, and Bm = 40 mT (solid line). 



at 1 ^ Kr < 1/6 as 
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with = For the s-wave bound state with m = 

this is the textbook asymptotic [H, [3] accounting for the 
conventional positive MR quadratic in small B . On the 
contrary, for orbitals with nonzero orbital momentum the 
wave function, Eq. (jlO[) . is linear in small B. 

If there is no time-reversal symmetry breaking, the 
states with the opposite direction of the orbital angu- 
lar momentum, m and —to, are degenerate, so that the 
linear term in the conductivity, a ~ am + (J-m cancels, 



aib) = a{0) 
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But even in this case the hopping conductivity, a{b) 
first increases with the magnetic field (negative quadratic 
MR) and only then decreases with B (positive quadratic 
MR), if Kp^/r < 2Am? — 4. Due to a large numeri- 
cal factor (24), this negative quadratic MR dominates 
in the whole region of realistic impurity densities for any 
nonzero to. Ions that carry a magnetic moment break the 
time-reversal symmetry and split m and —to states. Such 
zero-field splitting (ZFS) gives preference to the hopping 
via orbitals with a lower ionisation energy (positive m) 
providing the negative linear MR. In the extreme case of 
a ferromagnet with a frozen magnetisation, magnetoresis- 
tance to hopping via nonzero momentum orbitals should 
be highly anisotropic changing from linear and negative 
in the field applied parallel to the magnetisation to linear 
but positive in the opposite field, if the splitting due to 
exchange field is large enough compared with the tem- 
perature. 

At high electric fields the conduction in organic and 
inorganic insulators is often injection and/or ionization 
limited where carriers tunnel from extended states to a 
bound impurity level or visa versa under a potential bar- 



FIG. 3: (Color online) Transition from negative 
to positive OMAR (%) at room temperature in 
ITO I PEDOT / pentacene/Ca device at different bias 
voltages (symbols, Ref. [l3[) described with Ea. p4|l with 
r = 0.01, Bs = 6.7 mT, B,„ = 50 T (upper solid line), 
r = 0.0023, Bs = 3.5 mT, B„ = 10 T (middle solid line), 
and r = 0.01, Bs = 63 mT, = 5 T (lower solid line). 



rier shaped by the electric field [11[. The tunnelling ion- 
ization/injection rate, W can be described by a modified 
tunnelling ionization formula 12] fitting the numerically 
calculated ionization rates of atoms and ions over a large 
region of the electric field, F, 



Wm cx F- 
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where a and /3 are numerical constants, depending on a 
particular ion, and -F^ cx |eo — hujcjn/2\^/'^ is the charac- 
terstic ionization electric field depending on the magnetic 
filed in our case. Expanding Fi in powers of B yields 



Wm{F,b)^W{F,Q)e 



(13) 



where W[F, 0) cx F'^^ exp {-Fio/F - aF/F,o) is the ion- 
ization/injection rate without the magnetic field, FiQ cx 

Eq^^, and b = B / Biq is the reduced magnetic field with 
Bio = 4BQFFiQ/3{F^Q - aF^), which is a super-linear 
function of the electric field. Like in the case of the hop- 
ping conduction the injection/ionization conduction is a 
linear function of the weak magnetic field, resulting in the 
linear negative MR, if the orbital degeneracy is lifted. 

Finally, let us address puzzling experimental observa- 
tions of negative arid positive low-field MR in a number of 
organic materials [13l - [l5| . OMAR reaches 10% at fields 
on the order of only 10 mT, and can be either positive or 
negative, depending on operating conditions. The Zee- 
man energy does not account for the observed OMAR at 
ambient temperatures since it is too small, fisB sa 10 
mK, at a field of 10 mT. OMAR in hole-only devices 
[Tsl . [T^ rules out exciton-based mechanisms as an expla- 
nation. An alternative model involving spin-dependent 
bipolaron formation in deep potential wells has been pro- 
posed [13]. Nevertheless the origin of OMAR is still de- 
bated. Finding a convincing explanation of OMAR is 
important for understanding the basic transport mecha- 
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nism of organic insulators, which are used in molecular 
spintronics [TgI, [l7|. 

Here we propose a more general model describing var- 
ious classes of materials characterized by hopping trans- 
port regime, which accounts also for OMAR. We sug- 
gest that the hopping conductance could be a combina- 
tion of hoppings via conventional s-wave centers and via 
non-zero angular momentum orbitals. There is exper- 
imental evidence for paramagnetic centers and ZFS in 
polymers, in particular in Alq3 [T8| . More recently (su- 
per)paramagnetic susceptibility and ferromagnetic nan- 
oclusters have been reported in Alq3 [l^. The conven- 
tional hopping magnetoconductivity is described by the 
familiar exponential law as — aso exp[{—B /Bg)^], 
where Bs depends on the localisation radius and the 
density of s-wave centers, and aso is the zero- field con- 
ductivity (see also Ea. ((TT|) with m — 0). The conduc- 
tivity via non-zero angular momentum orbitals split in 
zero field is linear, Eqs. (|llll3p and FiglU so one can ap- 
proximate it as (T,„ = o'„io(l + B/Bjn), where amo and 
Bm = Bo/{2\m\Kp) or Bm = Bio/\m\ do not depend on 
the magnetic field. As the result the combined weak-field 
MR is described by the following simple expression. 



MR: 



-B/Brn + r[l^exp{-ByB^^)] 
1 + B/Brn + r cxp{-B^/Bf) ' 



(14) 



which can be readily compared with experimental data 
for sufficiently weak magnetic fields, B < Bm (here 
r = (Tso/cmo)- As one can see in Fig HI the theory de- 
scribes reasonably well the large negative OMAR mea- 
sured at room temperature in ITO / PEDOT / Alqi/Ca 
device at the bias voltage 14 V (and other voltages in 
Fig. 11 of Ref.flS]) using r,Bs,Bm as fitting parameters 
in the weak-field region, B < Bm- The low value of Bm 



points to a dominating role of hoppings via defects with 
rather shallow bound states (the binding energy on the 
order of a few K) in this device. Some empirical laws 
p^ . in particular cx — [BI{B-\- constant)]'^ also gives ac- 
curate agreement, so that the good fit might be coinci- 
dental. However, as noticed in Ref. _13:] simple fitting 
functions can fit the data only if one stays away from the 
transition region between negative and positive OMAR. 
Remarkably, as illustrated in Fig|3l simplified Ea. p4)) 
accounts for the cumbersome MR also in those organic 
structures, which show the transition from negative to 
positive MR. Compared with ITO / PEDOT / Alqi/Ca 
device, the ITO / PEDOT /pentacene/Ca device in Fig. 
14 of Ref.[l3] has deeper bound states. The electric-field 
behaviour of Bm and r are in agreement with the injec- 
tion/ionization conduction, Ea. (jl3p with a positive /3. 

In conclusion, we extended the conventional theory 
of hopping magnetoresistance to hoppings via non-zero 
orbital momentum states. The asymptotic 2D and 3D 
solutions of the Schrodinger equation show unusual lin- 
ear expansion/shrinking of the bound state with posi- 
tive/negative magnetic quantum numbers far away but 
not too far from the point defect in the magnetic 
field. Our theory accounts for an extraordinary nega- 
tive OMAR and for the transition to a more ordinary 
positive MR in disordered 7r-conjugated organic materi- 
als. Negative MR in some inorganic semiconductors may 
be also reanalyzed in the framework of the theory. 
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